Polaronic distortion and vacancy-induced magnetism in MgO 
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The electronic structure of the neutral and singly charged Mg vacancy in MgO is investigated 
using density functional theory. For both defects, semilocal exchange correlation functionals such 
as the local spin density approximation incorrectly predict a delocalized degenerate ground state. 
In contrast functionals that take strong correlation effects into account predict a localized solution, 
in agreement with spin resonance experiments. Our results, obtained with the HSE hybrid, atomic 
self-interaction corrected and LDA+C/ functionals, provide a number of constraints to the possibility 
of ferromagnetism in hole doped MgO. 



In recent years, the so called "d° magnets" have at- 
tracted great interest within the magnetism community 
and have been the topic of many theoretical and exper- 
imental investigations. Broadly speaking, the term "d° 
magnets" refers to a class of materials which, lacking 
any magnetic ions with open d or / shells, should in 
principle not be ferromagnetic, but nevertheless exhibit 
signatures of ferromagnetism often with a Curie temper- 
ature exceeding 300 K-i. Defective graphite 2 - and gold 
nanoparticlcs capped with organic molecules 3 arc clas- 
sical examples. However, the most widely studied class 
of d° magnets is that of un-doped or light-element doped 
medium to wide gap oxides. The first experimental study, 
which boosted research in this field, was on un-doped 
HfO^- Since then, there have been similar reports of fer- 
romagnetism for numerous other oxides such as TiGY 5 , 
SnOA CeO^, ZnQ£, CuO-A MgOl2«ii»i2 etc. Gener- 
ally, the ferromagnetism is observed in highly defective 
samples leading to the expectation that the magnetism 
must be somehow defect-related. Nanostructures are also 
often reported ferromagnetic and it has been argued that 
ferromagnetism can be a general property of nanopar- 
ticlesi^. Nevertheless the precise mechanism behind the 
reported high temperature ferromagnetism is not well un- 
derstood and is a subject of open debate. 

On the theoretical front, ab initio calculations based on 
Densisty Functional Theory (DFT) have been widely em- 
ployed to study d° magnetism in a wide range of systems. 
Following the early work of Elfimov et al on CaGi&, and 
in parallel with many experimental studies, a number of 
theoretical publications, proposing new d° ferromagnetic 
materials, have appeared over the past few years. In most 
of the d° oxides proposed, the formation of magnetic mo- 
ments and the resulting ferromagnetism is attributed to 
spin-polarized holes residing on cation p-orbitals either at 
vacany or impurity site a 16 ' 17 . Critically however, almost 
all such predictions are based on DFT calculations per- 
formed with local approximations of the exchange cor- 
relation potential (LSDA or GGA). These notoriously 
suffer from the spurious self-interaction, whose conse- 
quence is that of over-delocalizing the charge density. It 
is then not surprising that most of the calculations return 
a metallic (often half-metallic) ground state and usually, 
extremely large magnetic interaction. It has been shown 
previously that correcting for self-interaction can lead 



to the holes states localizing as a consequence of lattice 
distortions and the magnetic coupling is then reduced 
drastically or even completely suppresse d 18 ' 19 . This is 
consistent with the very well established notion that 
holes bound to acceptor defects form small polaron a 20 ' 21 . 
More recently, the electronic properties of cation vacan- 
cies in many Zn chalcogenides was investigated in detail 
by Chan et al^. It was shown that the metallic band 
structure predicted by LDA is changed into an insulat- 
ing one as soon as a hole-state potential operator—, which 
increases the splitting between occupied and unoccupied 
anion p-states, is applied. Pardo et al2* investigated the 
effect of on-site U corrections on the p-orbital magnetism 
in N doped alkaline earth monoxides and showed that 
for realistic values of the [/-parameter, the ferromagnetic 
coupling is drastically diminished. 

Without entering into the issues connected with per- 
colation 2 ^, already the fact that holes around acceptor 
defects form small bound polarons challenges the simple 
explanation of ferromagnetism due to carrier-mediated 
long-range coupling between magnetic moments residing 
on cation vacancies. The prototypical example of the 
cation vacancies in MgO clearly illustrates this point. Re- 
cently, numerous theoretical papers have proposed that 
Mg vacancies can induce room temperature ferromag- 
netism in MgOii^^. Based mainly on DFT LSDA and 
GGA results, Mg vacancy centres are reported to pro- 
duce half-metallic shallow acceptor states at the top of 
the valence band that lead to long range ferromagnetic 
coupling via a double exchange mechanism. 

However, the fact that a number of electron spin reso- 
nance (ESR) experiments dating back to the 1970s have 
clearly established the localized polaronic nature of Mg 
vacancies in Mg 26 i 27 seems to have been forgotten. As 
pointed out by Stoneham et al£^, the cation vacancy can 
be experimentally found in two charging states: neutral, 
V Mg, and singly charged, Vj^ g . ESR measurements es- 
tablish that V~M g s, in both the charging states, are deep 
traps exhibiting states in the gap as opposed to shallow 
levels at the top of the valence band. In the case of V Mg , 
the two holes from the vacancy localize completely on 
two adjacent oxygen sites. In contrast, in the case of 
V Mg , the only hole is completely localized on one single 
oxygen anion with an accompanying bond distortion. 

The availability of ESR experimental data, establish- 
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V Mg 


v ME 


LSDA 


2.19 (6) 


2.21 (6) 


GGA 


2.233 (6) 


2.24 (6) 


ASIC 


2.11 (2), 2.22 (4) 


2.13 (1), 2.24 (1), 2.25 (4) 


HSE 


2.111 (2), 2.253(4) 


2.106 (1), 2.225 (1), 2.253 (4) 



TABLE I: Calculated bond lengths (in A), for V Mg and V Mg 
as calculated by different DFT functionals. Bond lengths cor- 
respond to the distance from the center of the Mg vacancy site 
to the O ions in the co-ordination octahedron. The number 
in brackets is the number of bonds of a given length. 
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ing the nature of Vyig m MgO, sets a stringent test for 
DFT predictions obtained with different exchange cor- 
relation functionals. In particular the results of scmilo- 
cal LSDA/GGA can be compared to beyond LDA ap- 
proaches. In this work we look at three such methods. 
The first is the HSE hybrid-DFT functional 28 , which rep- 
resents a new and promising tool to study many prop- 
erties of solid state systems, in particular wide-gap ox- 
ides^. This was recently shown able to describe correctly 
the interplay between electronic properties and lattice 
distortion in oxides as complex as BaBiOj^. The second 
is the atomic self-interaction correction (ASIC) scheme, 
which has been employed previously to study a number 
of oxide and nitride wide-gap materials 31 . The third is 
the well known LDA+U method^. 

LSDA, ASIC and LSDA+J7 calculations were car- 
ried out using a development version of the SIESTA 
code^ 3 .. Norm-conserving pseudopotentials with the 
following reference atomic configurations were employed. 
Mg([Ne]3s 2 ), ([He]2s 2 ,2p 4 ). A double zeta basis set with 
additional polarization functions was used for both Mg 
and O. At the site of the vacany, extra basis functions 
corresponding to Mg atom were introduced to minimize 
basis set errors. Supercclls containing 96 atoms and 
192 atoms were used for single defect and multiple 
defect calculations respectively. Ionic coordinates were 
optimized using a conjugate gradients algorithm until 
all the forces are smaller than 0.04 eV/A. The GGA 
and HSE calculations were carried out using the VASP 
package^. A 54 atom supercell was used to study the 
V°M g center, while a cubic 64 atom supercell was used 
to simulate the singly charged Vj^ g center. A planewave 
energy cutoff of 520 eV and a 3x3x3 T-centered fc-point 
mesh were employed. For the HSE functional, the 
default value AEXX=0.25 is used for the Hartree-Fock 
exchange mixing parameter and for the screening of 
the Fock exchange HFSCREEN=0.2 is used. Geometry 
optimization was carried out using conjugate gradients. 

We start our analysis by presenting the LSDA results, 
with the GGA ones being rather similar. The VM g site 
in MgO presents a perfect octahedral symmetry which is 
preserved in LSDA, even upon lattice relaxation. This 
is the case for both V°M g and Vj^ g . The bond lengths 
resulting from this isotropic relaxation arc presented in 



FIG. 1: LSDA electronic structure of V Mg o- (Left) DOS of 
a MgO supercell containing either V°M g or V Mg . (Right) 
Charge density corresponding to the minority spin holes at 
the top of the valence band for V°M g - Note that the hole 
density is spread uniformly over the O octahedron. 



table [T] The associated hole density always spreads uni- 
formly over all of the six coordinating O anions with p- 
orbital lobes pointing towards the center of the vacancy 
(see Fig.[T]). In the case of the V°M g a spin triplet, S = 1, 
is obtained, while the spin singlet S = cannot be stabi- 
lized. In contrast V^ g presents a S =1/2 ground state. 
For both the charging configurations the density of states 
(DOS) is half metallic as shown in Fig.Q]and corresponds 
to that of a shallow acceptor impurity with hole states 
situated at the top of the valence band. This picture 
however does not conform with the experimental ESR 
data referenced earlier. 
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FIG. 2: ASIC electronic structure of V Mg o- (Left) DOS of a 
MgO supercell containing either V°M g or V Mg . (Top right) 
Charge density corresponding to the holes of V°M g - Note 
that the two holes are completely localized on two opposing 
O ions. (Bottom right) Charge density corresponding to the 
single hole of V°M g - The hole density is localized on one O 
ion. 

In contrast, a qualitatively different description of the 
ground state emerges from HSE and ASIC, with the two 
methods being consistent with each other. Firstly, lattice 
relaxation around the vacancy site is seen to break the 
perfect octahedral symmetry. In the case of the (V°Mg) 
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center, two virtual bonds from the center of the vacancy 
to two opposing O anions are found to shorten with re- 
spect to four other bonds forming a plane perpendicular 
to the two shorter bonds. In the case of the V^ g site, one 
of the two previously shorter bonds is seen to elongate 
while the four fold symmetry of the bonds in the plane 
perpendicular to these bonds is retained (see table UJ. 
The symmetry breaking around the vacany site is also 
accompanied by a strong localization of the hole density 
on the O ions nearest to the vacancy forming small po- 
larons (see Fig. [5]). It should be noted that the symmetry 
lowering around the vacancy site may not occur sponta- 
neously during the course of the geometry optimization, 
since the octahedrally symmetric solution is a local min- 
imum. Therefore, in order to access the lower symmetry 
ground state it might be necessary to initialize the geom- 
etry optimization in a slightly symmetry distorted ionic 
configuration around the site. 

HSE returns us an energy gain due to the polaronic 
distortion of 360 meV for V°M g and 225 meV for V^ Ig , 
and therefore it establishes that the distorted geometry 
is considerably more stable than the octahedrally sym- 
metric one regardless of the charging state. The ground 
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FIG. 3: DOS calculated with HSE for a MgO supercell con- 
taining either V°M g or V^ g . The top two panels show the 
DOS for the polaronic ground state geometry (labelled GS) 
of the defects. The bottom two panels show the DOS corre- 
sponding to local minima solutions (labelled LM) which pre- 
serve the octahedral symmetry around the vacancy sites. 



state DOS obtained within ASIC and HSE is shown in 
figures [2] and G3 As opposed to the shallow acceptor 
states observed in LSDA/GGA, both V°M g and V^ g ex- 
hibit a deep acceptor level in the gap. This is located 
at ~1 cV (~1.25 eV) above the valence band maximum 
(VBM) for HSE (ASIC). It is interesting to note that for 
the symmetrically relaxed local minimum solution within 
HSE (see lower panels of figure [3]), a shallow acceptor 
DOS similar to the LSDA solution is obtained. Thus the 
deep acceptor ground state is intimately associated to the 
polaronic bond distortion around the vacancy. 

In order to further substantiate the HSE/ ASIC re- 
sults we also carried out some test calculations within 
LSDA+J7 mainly looking at the evolution of the ground 
state as a function of the strength of the U paramc- 
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FIG. 4: (Left) LSDA+[/ energy difference between the LSDA 
and ASIC geometries for V°M g and V^ g as a function of U. 
Note the crossover between the two geometries at U ~ 4 eV, 
when the polaronic ASIC geometry becomes the most stable. 
(Right) Ground state DOS calculated with LSDA+t/ as a 
function of U. 



tcr. The U correction in this case is applied to the p- 
orbitals of the O sublatticc. For a given value of U, 
the total LSDA+L 7 energies obtained at the LSDA and 
ASIC geometries presented in table fl] are compared. The 
energy difference AE = £(G L sda) - -E(Gasic), where 
Glsda (Gasic) denote the LSDA (ASIC) geometry, is 
shown in figure |4j As the strength of U is increased the 
bond distorted ASIC geometry is seen to be stabilized 
with respect to the symmetric LSDA geometry and the 
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crossover is around U = 4 cV for both V m s 
Also the nature of the defect ground state changes from 
shallow acceptor to deep acceptor on either side of the 
crossover (Fig. 0]). In wide-gap oxides, the value of U 
for the 0-2p orbitals is estimated to be in the range of 
5-7 eV2i. At U=5 eV, for the V° Mg (V^g) center the 
polaronic Gasic geometry is lower in energy by a robust 
^390 meV (160 meV) with respect to the Glsda geom- 
etry. Furthermore the DOS for V°Mg (Fig. 2J exhibits a 
deep acceptor state located at roughly 1.3 eV above the 
VBM which agrees well with the HSE results. 

We now look at the magnetic properties associated to 
Mg vacancies. For the V°Mg center, the spins of the two 
holes may be aligned parallel or anti-parallel leading ci- 
ther to a 5 = 1 triplet or a 5 = singlet. As discussed 
above, LSDA/GGA predict a stable 5 = 1 configuration. 
However, within HSE/ASIC, the energies of the 5=1 
and the 5 = solutions are found to be very similar, sep- 
arated only by 0.75 meV (2 meV) in HSE (ASIC). We 
have also calculated the inter-defect magnetic coupling 
interaction between two (5=1) V°M g centers in a 192 
atom supercell. For LSDA two sites separated by 8.34 A 
are ferromagnetically coupled with a spin-flip energy of 
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33 meV (the energy required to reverse the spin of one of 
the two vacancies). This suggests the stabilization of fer- 
romagnctism via magnetic percolation between the holes, 
provided they are formed at a sufficiently high concen- 
tration. In contrast ASIC predicts for the same supercell 
configuration, a tiny interaction energy of 1 meV. Thus, 
we find that as soon as the correct polaronic ground state 
is achieved for the holes, their magnetic interaction vir- 
tually disappears. 

Before concluding, we add an additional remark on 
the similarity between V _ Mg and Li impurity substitut- 
ing Mg (LiMg). Indeed, both are single acceptors located 
at the cation sublattice. For LiMg as well, LSDA returns 
a half-metallic ground state, with the hole completely 
delocalized on all the O ions coordinating the impurity. 
ASIC instead localizes the hole on one O ion in a dis- 
torted polaronic geometry. This is qualitatively identical 
to what found for V~Mg (the Li-0 bond lengths are: 2.31 
A(l), 1.93 A(l), 2.15 1(4)). Although we did not find 
an experimental reference to confirm the ground state of 



LiMg in bulk MgO, its nature as a small bound polaron is 
well established in the case of surfaces and beyond-LSDA 
electronic structure methods have previously been used 
to study this defect in the context of catalysis^. 

In conclusion, we have shown that semilocal 
LSDA/GGA can lead to qualitative failures in their 
description of hole centers in MgO. Beyond-LDA ap- 
proaches that are either self interaction free or effectively 
correct for it, are able to reproduce the experimentally 
observed polaronic ground state. Significantly, the large 
inter-site ferromagnetic interaction predicted by LSDA 
is also shown to be an artifact. Thus, the observed high 
temperature ferromagnetic signals in MgO cannot be ex- 
plained by a simple model where magnetic interactions 
between hole centers, in their ground state, percolate 
through the sample. 

This work is supported by Science Foundation of Ire- 
land and by the EU (ATHENA project). Computational 
resources have been provided by the Trinity Center for 
High Performance Computing and by ICHEC. 
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